Highly-anisotropic elements for acoustic pentamode applications 
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Pentamode metamaterials are a class of acoustic metafluids that are characterized by a divergence 
free modified stress tensor. Such materials have an unconventional anisotropic stiffness and isotropic 
mass density, which allow themselves to mimic other fluid domains. Here we present a pentamode 
design formed by an oblique honeycomb lattice and producing customizable anisotropic properties. 
It is shown that anisotropy in the stiffness can exceed three orders of magnitude, and that it can be 
realistically tailored for transformation acoustic applications. 

PACS numbers: 43.40.+S, 81.05.Xjm, 62.30.+d 



Introduction. Metamaterials comprise a relatively new 
area of research, which utilizes a wave-functionalized mi- 
crostructure to yield unconventional effective material 
properties. With its origin in microwave optics [1], do- 
mains of study have since expanded to include near- 
infrared and visible optics radio-frequency electron- 
ics Q and magnetics with both linear and nonlinear 
phenomena being addressed. 

Many of the exotic phenomena observed in optical 
metamaterials, such as negative refraction Q and near- 
zero permittivity Q, have inspired analogous studies in 
acoustics. Both acoustic Q and elastodynamic [13] neg- 
ative refraction have been experimentally observed. Nar- 
rowband broad-angle negative refraction was achieved 
using simultaneous negative mass density (dipole) and 
stiffness (monopole) [ll-13|, or by novel anisotropy 
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utilizing perforated plates. 

Conversely, coordinate transformation methods 
have generated many new acoustic device concepts that 
do not require material negativity, such as lenses [l6| . 
beam sphtter s [iTII , black holes [18| and scattering reduc- 
ing cloaks 3" 22] ■ Such devices are defined as metafluids, 
which are effective materials with unconventional fluid- 
like properties whose particular bulk realization typically 
requires an anisotropic mass density, or in some limiting 
cases a tensor elasticity for flexural waves in thin plates 
[23} . Experimental demonstration has been sparse, with 
most studies relying on a superlattice approach of alter- 
nating isotropic layers 2^. However, such an approach 
is difficult to realize and limited by the so-called mass 
catastrophe, which requires infinite mass density in the 
effective material profile. 

An alternative approach is to generalize the conven- 
tional stress strain relationship to include pentamode 
elastic metamaterials. Pentamode materials 2^27] are 
metafluids that support five easy infinitesimal strains (i.e. 
there is only one non-zero eigenvalue of the elasticity ten- 
sor which is of a pure pressure type), and satisfies the 
invariance of the governing equations by virtue of main- 



taining a harmonic transformation. Pure pentamodes, in 
general, have an isotropic density and anisotropic stiff- 
ness with a negligible shear modulus. Recently, isotropic 
pentamode materials have become experimental reality, 
and it was shown that the structure's effective bulk mod- 
ulus exceeded the shear modulus by three orders of mag- 
nitude [28]. However, it has yet to be reported that 
anisotropic pentamode metamaterials can be realistically 
implemented for specific applications, since an elastic 
solid with a zero shear modulus would have no stabil- 
ity and immediately flow away. 

In this Letter, we show that an oblique honeycomb 
lattice can be utilized as a simple yet versatile build- 
ing block for pentamode device construction, which ex- 
erts highly anisotropic control over sound waves. The 
method presents a distinctly different approach to acous- 
tic metamaterials, in that it does not require the difficult 
to achieve high value anisotropy in the effective mass den- 
sity in addition to removing frequency bandwidth prob- 
lems associated with inertial metafluids. Potential ap- 
plications include extraordinary scattering reduction and 




FIG. 1. (Color online) Schematic of an anisotropic two- 
dimensional pentamode microstructure, with one unit cell 
highlighted. Bulk effective properties were tuned with varia- 
tion in cell geometry including joint diameter D, strut angle 
/? and inner-strut angle 7. The unit cell angle a is obtained 
by a = tan-^(l -I- tan(/3)), with -15° < /3 < 30°. ai is held 
fixed throughout. 
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FIG. 2. (Color online) Properties of an anisotropic pentamode 
metafluid. (a) Anisotropy factor ^ as a function of unit cell 
geometric parameters, using the constituent properties of steel 
and lattice orientation as in Fig. [1] Pentamode factor Ax in 
the stiff direction (b) and Ay weak direction (c). Note scale 
difference between (b) and (c), and that the y-axis in (b) 
and (c) is the same as in (a), (d) Shows the behavior of 
quasi-modes as the design becomes more anisotropic (small 
j3) at high propagation angles {6, which is measured from 
the horizontal axis). The volume fraction for the range of 
parameters covered in (a-d) spans from approximately 1.5% 
to 64%. 



arbitrary wave manipulation, low loss acoustic delay lines 
and phase controlled logic gates [30|. 

Anisotropy in pentamode metafluids. We consider elas- 
tic wave propagation in a microstructure having the gen- 
eral characteristics presented in Fig. [T] For simplicity 
we present our results in a two-dimensional (2D) plain 
strain space, however, the analysis can straightforwardly 
be extended to three dimensions. When the propagating 
wavelength is much larger than the lattice constants ai 
and a2, the microstructure can be assigned bulk homo- 
geneous eflFective properties. As shown in Fig.[T]the strut 
angle /? can be tuned to create isotropic (/? 30°), highly 
anisotropic (/3 ~ 0°) and re-entrant (/3 < 0°) materials. 
Hence, based upon a targeted property profile, a gradient 
structure can then be designed having smoothly varying 
properties related to the local changing microstructure. 



The condition that ai remain constant throughout does 
not effect the generality of the results presented. 

To demonstrate the properties of the structure, finite 
element analysis (COMSOL Multiphysics) was used to 
determine the structural eigenfrequenies of each geomet- 
ric permutation, which has Bloch-Floquet periodic con- 
ditions prescribed on the unit cell boundaries. This is 
analogous to computing the elastic wave band structure 
where the first two eigenfrequencies correspond to quasi- 
shear (qS) and quasi-longitudinal (qL) modes, respec- 
tively (see below for discussion on quasi-modes). These 
two eigenfrequencies were selected at small values of the 
wave number, near the origin where the dispersion curves 
are linear, to obtain long wavelength quasi-shear Vqs and 
quasi-longitudinal VqL phase speeds. An anisotropy fac- 
tor ^ is then defined as the ratio of effective bulk modulus 
in the principle directions, 
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Similarly, a pentamode factor A, which assess the con- 



tribution of shear, is defined as 



(2) 



^qL 



where the subscript represents the propagation direction. 
Inversion methods |3l| have also been widely used to re- 
trieve metamaterial effective properties, however such a 
method retrieves only a single modulus and is sensitive 
to termination conditions in the finite metamaterial slab 
used in the simulation. 

In the above equations it is assumed that the effective 
dynamic mass density p is isotropic and equal to the vol- 
ume averaged mass density (and hence is removed from 
Eqs. (PIS])). However, this assumption must be justified 
since p can be anisotropic \^2\ . Isotropy was verified by 
calculating the global reactant force and acceleration on 
the unit cell for a given harmonic displacement in the 
principle directions. It was concurrently determined that 
p was frequency independent within the long wavelength 
limit (see Supplemental Material). 

Figure [2{a) displays the calculated anisotropy factor, 
using the material properties of steel for the lattice and 
void for the interstitials, as a function of the geometric 
parameters /3, 7 and D, where D = D/|ai |. Exponentially 
increasing anisotropy is seen as f5 approaches 0°, with ^ 
exceeding 10'^ for small 7 and D. Figures 2(b)-(c) display 
the calculated pentamode factor, for the same geometric 
conditions as in Fig. 2(a). When /3 is small, A itself can 
be anisotropic; the contribution to shear modes is negligi- 
ble in the x direction (Fig. 2(b)) and substantial in the y 
direction (Fig. 2(c)). This occurs because as fi decreases 



the compression modulus (^ 



^qL 



) also decreases in the 



weak direction, yet the shear modulus {'^ '^qs) remains 
constant. 
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In isotropic media, the polarization direction of the 
propagating wave is either paraUel or normal to the wave 
vector direction. In anisotropic media, however, this is 
in general not the case and waves propagate as quasi- 
longitudinal (qL) and quasi-shear (qS) modes. The angle 
between the propagation direction = k/|k| and the qL 
polarization direction is expressed as 



c-gL = cos i(nfc • h^l)- 



(3) 



Since qL and qS modes are orthogonal (i.e. Ugs-tigL = 0) 
Ok-qS = Ok-qL + f , and for pure modes 9k-qL = 0. It 
is seen in Fig. 2(d) that under the current design there 
is a very strong dependence of Ok-qL on /3, and that for 
small /? there exist an extreme resistance of the wave to 
displace out of the stiff (x) direction. In summary, there 
are two phenomena which have the potential to affect ma- 
terial performance from the pentamode ideal, and both 
are born out of the relationship between the stiff (stretch- 
dominated) and weak (bending-dominated) directions of 
the lattice in the extreme anisotropic cases when /? is 
small. One (Fig. [UJc)) involves an increasing role of the 
shear mode in the weak {y) direction due to a decreasing 
compression modulus, and the other (Fig. [Sfd)) involves 
extreme deviation from pure propagating modes which 
tend to polarize along the stiff (x) direction for qL and 
weak {y) for qS. 

Device applications. To demonstrate the efficiency of 
our method to devise acoustic metafluids, we employ it 
to design two different transformation acoustic devices. 
The first one is based on the concept of an acoustic mi- 
rage, where an observer hears an echo from a distant 
wall, whereas in reality the echo originates from a much 
closer boundary. Norris 27| theoretically showed that a 
2D mirage can be implemented using a pure pentamode 
metamaterial having orthotropic symmetry in the stiff- 
ness tensor. This metafluid, of length a, mimics a larger 
fluid domain, of length a + d, and is both impedance 
matched and has equivalent acoustic travel times, such 
that to an external observer the two domains are indis- 
tinguishable (see Fig. S}. The mirage density and bulk 
modulus are expressed as (2?} 











(4) 



where po and i^o ai'e the density and bulk modulus of 
the background fluid, in the present case water. Equation 
(j4|) describes a fluid with an orthotropically symmetric K, 
and the mirage device can now be designed by choosing a 
constituent material and adjusting the unit cell geometry 
to give the required properties. 

Figure [Sja) shows the band structure of a unit cell (see 
inset of Fig. El^a)) which approximately satisfies Eq. ^ 
for a mirage parameter of a = d, as a function of the 
normalized frequency Q. = 27r/|ai|/co, where / is the cir- 
cular frequency and cq is a reference sound speed. For 
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FIG. 3. (Color online) (a) Normalized band structure along 
the principal directions of propagation for the unit cell used 
in the acoustic mirage. Inset shows the unit cell and first 
Brillouin zone with the irreducible part shaded in gray, (b)- 
(c) Equifrequency contours for qS and qL modes at selected 
frequencies (contours not calculated to zone boundaries). 
The dotted curve depicts the wave vector k direction, and 
anisotropic behavior is depicted by the difference between k 
and group velocity vector Vg 



this choice of geometric parameters, Kn = 0.5052Ko, 
K22 = 2.OK0 and p = 2.037/90, as can be seen in the band 
structure. With using silver as the constituent material, 
^ = 0.24, = 0.96 and Aj, = 0.99 (note that the unit 
cell of Fig.|3]is rotated 90° as compared to that in Fig.[T]). 
The device has a sufficient amount of shear for structural 
stability yet still retains pentamode characteristics. Ad- 
ditionally, a wideband shear-wave directional band gap 
is noticed for propagation in the weak (x) direction. Fig- 
ure[3l^b-c) shows the equifrequency contours of the qS and 
qL modes at selected frequencies. At £7 < 0.01 the qS 
modes are approximately isotropic, with deviation from 
circularity of less than 5%. For higher Q, the qS modes 
have both significant increasing anisotropy and disper- 
sion, as indicated by the relationship between group ve- 
locity and wave vectors inside the lattice. The equifre- 
quency contours for the qL mode, conversely, maintains 
the required ellipsoidal anisotropy over a broad range of 

n. 

Using a unit cell with the specific properties obtained 
above, an acoustic mirage device is constructed as de- 
picted in Fig. m Figures ID^a) and IDJc) respectively com- 
pares the scattered acoustic fields for the mirage and the 
non-mirage case (the larger domain, termed the control). 
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FIG. 4. (Color online) Acoustic mirage device for a mirage parameter of a — d, showing the real part of the scattered pressure 
field, using water as the background fluid. Shown is the case for frequency = O.I and an incident angle Oi = 43.4°, geometric 
parameters are those from Fig. [S] An observer at a; = sees the same acoustic field (amplitude and phase) with the pentamode 
metafluid (a) as without (c). For reference, a domain (b) is shown with the same length as the mirage. A plane wave travels from 
left to right, with imposed zero-reflection condition on left boundaries and zero displacement condition on right boundaries. 
Top and bottom boundaries have Bloch-Floquet conditions. 



and they are seen to have excellent agreement. A refer- 
ence domain that is the same size as the mirage domain 
is also shown in Fig. |3^b). Since the pentamode design 
is not based on narrow resonances, the mirage device 
will also operate over both a very wide frequency band- 
width and all angles of incidence. Figure [SJa) and \5lh) 
summarizes the difference of the scattered pressure field 
amplitude dp and phase Scj) between the mirage and con- 
trol, with excellent agreement seen. This is compared to 
Figs. EJc) and[5jd), which depend strongly on the inci- 
dent angle and frequency between the two domains. 

The above outlined design method also has the flexi- 
bility to realize more complicated material profiles. To 
demonstrate this, we use the method to realize the es- 
timated material property values for a cylindrical pen- 
tamode scattering reduction layer for use in an aqueous 
environment. A graded anisotropic honeycomb lattice is 



used to simulate the transformed space, which reduces 
the actual size of the cylinder's radius a to virtual radius 
6. One possible mapping function (though not unique) 
describing the transformation layer is ^33] 



fir) = (62-a2)-i ib^-a6)r 



r 



(5) 



where r and h are the radial coordinate and outer radius 
of the transformation layer, respectively. Equation ([S]) 
maintains both radial impedance and orthogonal wave 
speed matching with the external fluid. Figure [6] com- 
pares the ideal pentamode material properties and val- 
ues obtained from the composite metafluid for the case 
oi 5 ~ a/Z and h = 2a, using steel as the base material. 
In 2D, the pure pentamode properties can be obtained 
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FIG. 5. (Color online) Difference of scattered pressure am- 
plitude (5p) and phase (50) of the mirage (a) and (b), and 
reference domain (c) and (d), evaluated at a point location 
far from the water-metafiuid interface, a; = 0. Plotted as a 
function of normalized frequency SI and incident angle di. 



FIG. 6. (Color online) Designing material properties for 
transformation acoustics scattering reduction. Solid lines de- 
note ideal pentamode properties and symbols are properties 
achieved under the current design. All values are normalized 
to water. 
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through [ii] 





fir) r 






-p{r) 


r 



where the over hne represent normahzation to water val- 
ues. Excellent approximation is demonstrated. For all 
cases of the designed properties, pentamode behavior is 
achieved, = - 1. At r = a, in Fig.lH i?„ = 0.16, 
K±^ — 6.33, and p = 0.70. An inertial metafluid ver- 
sion of the transformation layer using Eq. ([5]) produces 
(at r = a) 7? = 1.42, = 0.157, and p„ = 6.33. Us- 
ing a superlattice approach to realize these inertial val- 
ues requires alternating isotropic layers with properties 
Pi = 12.58, 'Ki = 6.29 and -p^ = 0.079, K2 = 0.80. Real- 
izing such a composite material is prohibitively difficult. 

Finally, it is noted that the approximation of the ideal 
properties in Eq. ([6]) are for a cylindrical coordinate sys- 
tem, whereas the unit cell analysis is done in a recti- 
linear system. However, it is the point of this work to 
demonstrate the feasibility of the design to approximate 
unique property values required in pentamode applica- 
tions, the design of a gradient curved structure based 
upon an anisotropic honeycomb lattice is beyond the 
scope of the current work, and will be presented in a fu- 
ture publication. Furthermore, since the current method 
is fully scalable, as the size of the unit cell is decreased 
the local difference between the curvilinear and rectilin- 
ear coordinate systems will be reduced. 

Conclusions. A versatile method to design and charac- 
terize highly anisotropic pentamode elements for trans- 
formation acoustics was presented. The design method is 
based upon customizing the geometric parameters of an 
oblique honeycomb lattice to target particular transfor- 
mation acoustics applications. In the low frequency limit, 
the metallic structures act as acoustic metafluids with 
pentamode properties, and have negligible shear mod- 
ulus relative to its anisotropic fluid-like stiffness. The 
simplicity of the design (easy to fabricate) coupled with 
the effectiveness of performance shows promise to realize 
concepts hitherto constrained by inertial-based metama- 
terial constructions. 

This work was supported by the Office of Naval Re- 
search. 
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SUPPLEMENTAL MATERIAL 



To verify isotropic mass density, the global reactant 
force Fi and acceleration Ui is obtained by averaging 
the local surface traction Tj — niUij and acceleration 
iii on the unit cell's external boundary E, where Ui 
and Uij are the local surface unit normal vector and 
stress tensor, respectively. These can be expressed as. 



Fj = ^ J^maijdr and Ui = j J^iii dr, where A,£ and 
dr are the area of the unit cell, total path length of in- 
tegration and differential length of integration along E, 
respectively. The effective dynamic anisotropic mass den- 
sity is therefore obtained simply as pi = Fi/ai. Figure [7] 
compares p and pi {i — x,y), using steel as the base ma- 
terial, and it is seen that they are approximately equiv- 
alent. Pi is also frequency independent within the long 
wavelength limit. 
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FIG. 7. (Color online) Comparison between volume averaged mass density p and anisotropic mass densities px and py. 
Equivalence is shown within the long wavelength approximation. 



